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Abstract 

In this paper, we restudy the Green function expressions of field equations. 
We derive the explicit form of the Green functions for the Klein-Gordon equation 
and Dirac equation, and then estimate the decay rate of the solution to the linear 
equations. The main motivation of this paper is to show that: (1). The formal 
solutions of field equations expressed by Green function can be elevated as a 
postulate for unified field theory. (2). The inescapable decay of the solution of 
linear equations implies that the whole theory of the matter world should include 
nonlinear interaction. 
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1 Introduction 

Since the Green function or propagator of the wave equation was introduced to calculate 
the scattering matrix by Stuckelberg[l] and Feynman[2], it has become an important 
tool of quantum field theory [3, 4]. Green function is also a useful tool for asymptotic 
analysis and L°° estimation of fields. The Feynman propagators are in the form of 
Fourier integrals with the fields being approximately treated as plane wave. In some 
old textbooks[5, 6], the explicit Green functions of the wave equations were expressed 
by Bessel function. In this pedagogical paper, we restudy this problem in detail. We 
derive the Green functions of the Klein-Gordon equation and Dirac equation as normal 
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functions, and express the general solutions of the field equations in integral form, then 
use this form to estimate the decay rate of the linear fields. 

The main motivation of this paper is to show that, (1). the solutions of field equation 
expressed by Green function are more naturally and visually display the correlation of 
fields than differential equations, and the expression of integral form are more flexible 
and comprehensive. This feature of the expression enables us to elevate it as a postulate 
for unified field theory. (2). The estimation of decay rate of linear fields shows that, the 
linear field theory is not complete enough to describe the stability of the fundamental 
particles, A whole theory of the matter world should include nonlinear interaction. 

2 Preparation 

In this paper we only consider the distribution defined in R + = [0, +oo). By straight- 
forward calculation, we have 

Lemma 1. For p G R + , in the sense of distribution we have 




(2.1) 



where H(p) is the Heaviside function, and 5(p) is the Dirac-5 




1, p > 
0, p = 




(2.2) 




dp 
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where J n {p) {n — 0, 1, 2, • ■ •) is the Bessel function 
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Lemma 3. Suppose m > 0. Let w = \/m 2 + p 2 and 




(2.5) 
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then in the sense of distribution, we have 



urn 2 



K l (t,r) = —iJ (p)6(p)-J 1 (p)H(p)}, p = mVW^ (2.6) 
zp 



(2.7) 



Proof. Let 

w = mcoshz, p = msmhz, 
t = — cosh P, r — 2- sinh P , 
then we have 

K i = 2 P skihg S-oo thz ■ |[cos(wt - pr) - cos(wt + pr)]dp 

= I^K^ I-oo sinh z l cos (P cos H z ~ 0) -cos(pcosh(z + P,))]dz 

= 4^h|/^oCOs(pcosh S )[sinh( S + e)-sinh( S -e)]^ (2-8) 

= — J °° cos(p cosh s) cosh sds 

= ff?COs( P T) 7 ^= l dT=^I>(p), 

where I is defined by (2.3). By J' Q = —J±, we prove the lemma. 
The following lemmas can be checked directly. 

Lemma 4. For all real vector p — (pi,P2,P3), denote p = Jp\ +p| +P§ an d 

-i / pi-ip2 Pi+ip2 \ 

B 1 (p) = -=( fl!L_ /^H), (2.9) 

then Bi is a unitary matrix satisfiing 

B^ = B+, B^a-pB 1= pa 3 , (2.10) 

where a are Pauli's matrices. 

Lemma 5. B\ is defined by (2.9), denote 

Q = m'y + a-p, B 2 = diag(B l , B x ), (2.11) 

7= (o -/)• s= (l o)' (2A2) 

then we have 

B 2 1 QB 2 = m 1 + paz. (2.13) 



Lemma 6. Let w = \/ m 2 + p 2 , 



( 



B, 



'Iw 



y/w—m 









y/w—m 

yjw — m 
V y/w — m 



y/w+m 





p 



y/w+m 





a/ic + m 

y/w + m J 



(2.14) 



3 



then we have 

B- 1 (m^+pa 3 )B 3 = wy. (2.15) 

Lemma 7. Let Q = B2B3, then 

Q = wQ 1 Q~ 1 or Q~ 1 eQ = w-f J (2.16) 
where all parameters are defined as above. 

3 The Green functions of field equations 

(Fl). For the initial value problem of the Klein-Gordon equation 

d 2 u = Au — m 2 u, with u(0,x) and ut(0,x), (3-1) 

the formal solution represented by the Green function is given by 

u(t,x) = I \K(t,r)u(0,y) + G(t,r)u t (0,y)}d 3 y, (3.2) 

where r — x — y. 

Lemma 8. In the sense of distribution, for (3.2) we have 

G(t,r) = — — / — sin wt sin pr dp, (3.3) 
27T 2 r Jo w 

I roc 

K(t,r) = d t G = — — / p cos wt sin pr dp. (3.4) 
zn^r Jo 

Proof. Denote u(0,x) = (p(x), u t (0,x) = ip(x), Make Fourier transformation to 
(3.1), we get 



where 



= — {p 2 + m 2 )u = —w 2 u, 
u(0,p) = <j>(p), u'(0,p) = <0(p), 



u(t,p) =±^J R3 u(t,y)e-^ i d?y, 
u(t, %) = 7y=f Xrs u(t,p)e 3 -&d 3 p, 



(3.5) 



(3.6) 



and so on. The solution of (3.5) is given by 



u — cos wt<p -\ smwtih. (3.7) 

w 



Substituting (3.7) into (3.6), we get 



Gf(t,r) = — L-/ - sin wte^-^d 3 p, (3 A 

v ' ' (2tt) 3 Jb? w F} k 
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K(t, r) = -—-^ cos wte^-^d^p = d t G, (3.9) 

Let f-p — pr cos 9, integrating (3.8) and (3.9) in the spherical coordinate system proves 
the lemma. 

By Lemma 3 and (3.3), we get 

Lemma 9. Let p = m\Jt 2 — r 2 , then we have 

777 

G(t,7) = —j=== [J (p)6(p) - Mp)H(p)}, (3.10) 



By Lemma 9, we have 

Theorem 1. The solution of (3.1) can be expressed by the following integrals 

u(t, x) = u (t, x) — Ui(t, x), (3-11) 
where uq is the solution for the case m = 0, 

"»= 9 '(4sL«) + isL«- (3 - 12) 

S y stands for the surface r = t, 



Proof. For simple, we set = 0. Then by Lemma 8 and 9, we have 

1 Jo* 7P=?^(p)r 2 dr £ sin Odd f^ ^ - Ul (t, x) 
± ft rd{-my/W^)J Q 8{p) £ sin 0d6 ifjdcp - Ul (t, x) 



4tt JO ^/t 
J_ 
An 

Jq sin QdO J 2?r ip(r = t)d<p — ui(t, x) 



In 



(3.14) 



The proof is finished. 

(F2). For the Dirac equation 

dti/j = —a ■ Vif) — mjiifj, with -0(0, x) = 4>(x), (3.15) 

by Fourier transformation we get 

Lemma 10. The formal solution of (3.15) is given by 

iP(t,x)= f A(t,x-y)<P(y)d 3 y, (3.16) 
Jr. 3 

where 

2777.7 1 ^1 , , 

A(t, r) = -^ Kl + —d t K t + —a ■ VK U (3.17) 
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and V = (d yi ,dy 2 ,dy 3 ), Ki{t,r) is defined by (2.5) or (2.6). 
Proof. By Fourier transformation we will get 

1 



A(t,r) 



where G = rwy + a ■ p. By Lemma 5 and 7, we have 



-en 



= cos wt — i sin wtQjQ -1 

= cos wt — ism wt{m^ + a ■ p), 

Substituting it into (3.18) we have 

A(t, r) = — - / [cos wt — i sin wt(mj + a ■ p)]e r '^ l d 3 p. 



Denote 



T 



f = r(sin 6\ cos (fx, sin 9\ sin <pi, cos 9i), 
/ cos 61 cos <pi cos 9i sin tpi — sin 9i \ 

— sin 9\ sin ipi sin 9\ cos <pi 
\ sin 9\ cos ipi sin 9\ sin ipi cos 9\ / 



then 



f = r(0,0, 1)T. 

Making transformation 

p = kT, k = fc(sin 9 cos ip, sin 9 sin ip, cos 9), 
substituting (3.23) and (3.24) into (3.20) we have 

a ■ r 

a- p— {(3i cos ip + j3 2 cos ip) k sin 9 H k cos 9, 

r 

where (3\, ft? are constant matrices, r-p — rk cos 9, 

A(t,r) = ^ k 2 dk £ sin 9[cos wt - ^ sin wt(mj + ^k cos 9)}e ikrcose d9. 

Integrating (3.26) with respect to 9 gives (3.17). 
By Lemma 3 and 10, we have 

Theorem 2. The solution of (3.15) can be expressed by 



where 



^0 = d t 



\4lTt Jr 



r=t 



4){y)dS l 



1 

Ant Jr 



r=t 



(mji + a ■ V)(f)(y)dSy, 



ipi = d t 



™ [ Ji{p)4>{y) 



An 



■<t y/t 2 — r 2 



rn 
An 



r<t 



Vt 2 



1(77172 + a ■ W)4>(y)d?y. 
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4 The decay rate of the linear fields 



Now we use (3.11) and (3.27) to estimate the decay of the solutions. In [7, 8], for wave 
equation we have 

Theorem 3. Suppose <f> £ W n+3,1 (R 3 ), tp £ W n+2,1 (R 3 ), then the solution of 



(4.1) 



(<9f - A)u = 0, 
u(Q, x) = 0, d t u(0, x) = if), 

has the following decay estimate 

C 

IK*, -)\\w^{R?) < YT^(H^llw i+3 ' 1 (« 3 ) + I Iv^l Iv^"+2.i(^3)), (n > 0), 
C 

\\d t u{t, •) | 1^-1,00(^3) < J— (||0|| W «+3,l (fl 3) + 11-011^+2,1(^3)), (n > 1), 

where C always stands for a positive constant independent of (/), if). 

For the Klein-Gordon equation (3.1), we have the following estimation. 

Theorem 4. Suppose £ W 3 ' 1 (R 3 ), ip £ W 2 ' 1 (R 3 ), then the solution of (3.1) has 
the following decay estimate 



(4.2) 
(4.3) 



VT+t 

Proof. By theorem 3, we have 

C 

MM) I < Y^~ t 

where Uq is defined by (3.12). Let 



C 

u(t,x)\ < : (||0||w3,i (fl 3) + \\lf)\\w^(B 3 )] 



W ^(R3) + \\l/)\\ W 2,l^), 



Mp) 



r<t 



Vt^ 



v(t, x) 

Since \J\\ < 1, for case t < 1, we have 

KM) I = \ti^drf Q d{lf r °°d s iPds 

where dQ = sin 6 d9 dtp. For t > 1, we have 



(4.4) 



(4.5) 



(4.6) 



(4.7) 



|f(t,f)| 



r 2 J\dr 



< 



rt dr 



t r 2 J\dr 



f n dQ r d s if)ds 



* 2 v /( i +r)(i-r) 

Combine it with (4.7), we get 



(4.{ 



C 



VTTt 
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(4.9) 



Similarly, we have 



0, 



Up) 



r<t 



Vt 



4>(y)d 3 y 



< 



C 



VT+t 



I^IIVK 2 . 1 ^ 3 )- 



(4.10) 



The theorem holds by (4.5), (4.9) and (4.10). 

Considering linear the Dirac equation also satisfies the Klein-Gordon equation (3.1 
so we have 

Corollary 5. For Dirac equation (3.15), we have the following estimation 

C 



t,x) < 



(4.11) 



5 Discussion and conclusion 

From Theorem 2 we learn that, the solution of a physical field system of the first 
order equation u(t,x) = (u±, ui ■ ■ ■ u n ) T can be generally expressed by the following 
functional 

G(u(x,t)) = [ dr f (F + K-V*u)d£, (5.1) 

Jt J(1(t) 

where G(u) is an invertible and differentiable function group of Uk(x, t), to is any given 
time less then t, U T6 [ i0ji )f2(r) is the dependent domain of u(x,t), 

F = F(u(£, r), r),C, r, x, t, t ), K = K{u{t t),v{& t, x, t, to), (5.2) 

all the elements of the group F and the vector matrices K are normal functions of 
u and v. v denotes the interaction of the environment. The physical meaning of 
(5.1) is that, the present state of the system is simply determined by its history and 
interaction. Combined with the principle of limited light velocity, Eq.(5.1) inversely 
implies the first order field equation 

A»d,u = f(u,v). (5.3) 

The reason to choose (5.1) instead of (5.3) as postulate to describe the field system is 
that, (5.1) is more flexible and comprehensive, and its visual physical meaning can be 
understood by a farm folk[9, 10]. 

By (4.11) we learn that, the solution of a linear wave equation will inescapably 
decay to zero. To explain the stability of the fundamental particles, a complete theory 
should include nonlinear Dirac equations [11, 12]. 
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